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Modes  in  chiral  core  planar  waveguides:  transition  from 
linear  to  circular  polarization 
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Abstract 

The  general  solution  for  modes  in  an  asymmetric  planar  waveguide  with  an  isotropic  chiral  core  is  given  in 
terms  of  a  pair  of  parameters  related  to  the  eccentricity  of  the  polarization  ellipse  for  the  transverse  electric  field. 
This  formulation  provides  insight  into  the  transition,  with  increasing  chirality  of  the  core,  from  TE/TM  modes  to 
right-handed  (RHC)  and  left-handed  (LHC)  circular  polarization  modes.  The  mode  properties  as  a  function  of 
waveguide  thickness  and  of  frequency  are  discussed  in  detail.  In  addition  to  the  usual  cutoff  at  a  minimum 
thickness  (and  maximum  wavelength),  left-handed  elliptical  modes  are  found  to  experience  a  cutoff  at  a  mode- 
dependent  maximum  thickness  (and  minimum  wavelength).  The  limiting  case  of  a  symmetric  waveguide  is  also 
discussed. 


1.  INTRODUCTION 

The  effect  of  chiral  media  on  electromagnetic  waves  has  been  the 
subject  of  investigation  for  two  decades,  with  a  good  deal  of  the 
initial  interest  centered  about  radar  wavelengths.1,2,3  In  keeping 
with  the  interest  in  that  part  of  the  electromagnetic  spectrum, 
early  studies  of  waveguides  containing  chiral  media  were  on 
dielectric  waveguides  with  conducting  boundaries.4,5,6’7  At 
optical  wavelengths  in  bulk  material,  one  of  the  important 
physical  effects  of  chirality  is  optical  activity.8  Chiral  polymers9 
are  potentially  excellent  materials  for  chiral  optical  waveguides 
provided  sufficient10  chirality  can  be  realized  in  a  polymer  with 
low  optical  loss  and  good  film-forming  properties.  Not 
surprisingly,  chirality,  or  handedness,  of  the  medium  in  all¬ 
dielectric  waveguides  has  been  found  to  introduce  considerable 
additional  mathematical  complexity  resulting  in  mode  properties 
that  can  substantially  differ  from  the  achiral  dielectric  case. 
Much  of  the  work  on  planar  chiral  dielectric  waveguides  in  the 
literature  deals  with  symmetric  boundaries, 11,12,13,14,15  in  which 
the  top  and  bottom  claddings  have  the  same  dielectric  properties. 
For  asymmetric  planar  waveguides,  perturbation  techniques  have 
been  employed1  16,17  or  approximations  have  been  made18  to 
simplify  the  mathematics. 

Here  we  develop  the  general  solution  to  the  asymmetric 
planar  dielectric  waveguide  with  an  isotropic  chiral  core  without 
making  simplifying  approximations.  The  allowed  modes  are 
elliptically  polarized  and,  in  our  formulation,  the  modal 
equations  contain  a  pair  of  parameters  that  determine  the 
eccentricity  of  the  polarization  ellipse  for  the  transverse  electric 
field.  The  solution  for  a  symmetric  waveguide  is  presented  as  the 
limiting  case  when  the  top  and  bottom  cladding  indices  are  equal. 
When  the  values  of  the  parameters  are  appropriate  for  transverse 
electric  (TE)  and  transverse  magnetic  (TM)  modes  in  the  limit  of 
vanishing  chirality,  the  modal  equations  are  easily  seen  to  reduce 
to  the  well-known  equations  for  the  TE  and  TM  modes  of  an 
achiral  waveguide. 

2.  FIELDS  FOR  CHIRAL  MEDIA 

We  use  exp(ico/)  time  dependence  for  Maxwell’s  equations, 
which  results  in  the  equations  involving  time  derivatives  taking 
the  forms 


VxE  +  /<wB  =  0  ,  V x H  -  icon  =  0  ,  (1) 

and  adopt  constitutive  equations  including  chirality  in  the  Drude- 
Bom-Federov  form,19 

D  =  e(E  +  y  VxE) ,  B  =  p(H  +  y  VxH) ,  (2) 

where  e  and  p  are  the  usual  permittivity  and  permeability  and  y  is 
the  chirality  parameter,  which  has  units  of  length.  The  wave 
equation  for  E  resulting  from  Eqs.(l)  and  (2)  contains  a  VxE 
term.  Wave  equations  in  chiral  media  that  reduce12,18  to  the 
Helmholtz  equation  can  be  obtained  by  employing  Bohren’s 
decomposition20  of  E  and  H  ,  which  may  be  taken  to  have  the 
form18 

E  =  — (F+  +  F~)  ,  H  =  i7J-(F+-F-)  ,  (3) 

2  2i\n 

with 

F*  s  E  ±  i-JJiJe  H  .  (4) 

From  Eqs.(l)  and  (2)  it  follows  that  the  vector  field  functions 
F1  satisfy 

VxF*=TV±F±  >  (5) 

where 

"±=7±s  ’  "«  ’  8  =k°n«r  ’ and  *°  =T" 

(6) 

All  example  calculations  in  the  figures  in  this  paper  are  done  for  a 
wavelength  of  633  nm.  Note  that,  in  terms  of  these  quantities,  the 
rotatory  power  in  bulk  material  is  p  ~k^ng8  for  8«1.  The 

vector  fields  F+  and  F-  separately  obey  the  Helmholtz  equation, 
V2F*  +  (A0«±)2F±  =  0.  (7) 
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Figure  1 .  Coordinate  system  for  an  asymmetric 
planar  waveguide  with  an  isotropic  chiral  core. 

In  a  bulk  medium,  F+  and  F"  are  plane  waves  with  right- 
handed  circular  (RHC)  and  left-handed  circular  (LHC) 
polarization,  respectively. 

3.  THE  ASYMMETRIC  PLANAR  WAVEGUIDE 

The  structure  of  the  planar  waveguide  under  consideration  is 
shown  in  Fig.  1 .  The  top  and  bottom  claddings  are  not  chiral  and 
have  refractive  indices  n0  and  ns ,  respectively,  while  the  core  is 

isotropic  and  has  permittivity  e  =  £qh|  and  chirality  Y-  The 

bottom  cladding  will  frequently  be  referred  to  as  the  substrate. 
The  waveguide  has  thickness  d  and  is  assumed  to  be  infinite  in 
the  ±x  direction.  The  notation  here  mostly  follows  that  of  Ref.  18 
with  a  notable  exception  that  here  we  use  the  symbol  y  for  the 
chirality  parameter  instead  of  p  used  there  in  order  to  clearly 
distinguish  between  the  chirality  parameter  in  the  constitutive 
equations  from  the  usual  symbol  p  for  rotatory  power. 

Waveguide  modes  propagating  in  the  +z  direction  have 

the  form 

F*  (y,  z)  =  V  ±(y) exp (~ik0neff  z) ,  (8) 

where  neS  is  the  effective  index  of  a  mode  with  propagation 
constant  /3  =  k§tieff  .  It  follows21,18  from  Eqs.(5)  that  the  y-  and 

z-components  of  *  can  be  found  from  the  jc-component,  and 

consequently  ¥  *  can  be  written 


nq 

^±(y)  =  ±1  dV~(y) 
k0nq  dy 


(9) 


where  y/~(y)  satisfy  the  Helmholtz  equation,  . 


—^-3—  +  k02(nq2  -n2ff)v±(y)  =  0. 


(10) 


Here,  q  -  ±  in  the  core  layer,  and  q  =  0,s  in  the  top  and  bottom 
claddings,  respectively. 

The  solution  to  Eqs.(10)  can  be  expressed  in  the  form 


v±(y)= 


A~  exp(-vy) ,  y  >  0 
fi±cos(«±y+^±) ,  0 <y<-d 
C±exp(w(y+^))  ,  y<-d 


(11) 


where 

«±s*0-Jn±~neff  ’  v-k0ijneff  ~n0  >  w=k0^neff  ~ns  ■ 

(12) 


A  mode  in  the  chiral  core  waveguide  consists6  of  RHC-  and 
LHC-polarized  waves  traveling  down  the  guide  with  the  same 
effective  index  neffy  but  total  internal  reflection  off  the  interfaces 
is  at  different  angles  of  incidence  0*  due  to  differing  refractive 

indices  n*,  i.e.y  sin 0*  - neff  jn+  and  cos#*  =  urj{k^n±) . 

There  are  four  such  circularly  polarized  waves  making  up  a  mode 
in  the  core:  one  pair  of  RHC-  and  LHC-polarized  waves 

exp( -/0 1 ) (s ± i p f ) exp( -/K f  -r)  (14) 

traveling  toward  the  top  interface  with  propagation  vectors 
Ki*  =  *0"±(°’  cos  0  * ,  sin  0  * ) ,  and  another  pair 


02=  exp(/0±)(s±ip2  )exp(-/K2  ■?)  (15) 

traveling  toward  the  bottom  interface  with  propagation 
vectorsK2±  =k0n±(Q ,  -cos 0*,  sin 0*).  The  caret  in  Eqs.(14) 

and  (15)  denotes  a  unit  vector  and  s  =  x ,  p^*  =Klt2:t  xx  . 

The  vector  functions  F*  are  related  to  these  circularly  polarized 
waves  by 

F*  +®f).  (16) 


The  boundary  conditions  requiring  continuity  of  the 
tangential  components  of  E  and  H  at  the  top  and  bottom 
interfaces  give  rise  to  8  equations,  which  for  completeness  are 
given  in  the  Appendix.  Using  trig  identities,  Eqs.(A9)  to 
(A  12)  can  be  written  as  linear  equations  for  B^in  ty* ,  B^os  in 
the  form 


M,,  (neff )  = 


B+  sin0  + 

(rieff) 

B+  cos<j>+ 
B~  sin#“ 
B~  cos <t>~  j 

=  0  , 

-1 

*0 

-l 

r0Gr0 

-1  - 

'0<T0 

i 

C+ 

-S  + 

C~ 

-5“ 

Cr+ 

-5r+  - 

cr 

Sr~ 

(17) 


where 


2 


(24) 


C~  =<r*cos(tt:trf)  +  sin(K±</)  , 
Cr~  -rsaj  cos(w~</)  +  sin(K“</) , 
S~  =af  sm(u~d)-cos(u±d)  , 
Sr±  =  rsaf  sin( urd)  -  cos iu^d) , 


X  X 

<t±=(1±5)—  ,  **=(115)—, 


„2 

_fo 


2  >and  r*“  2  • 


(18) 


001(0*  )  =  CTq 


*0*8 

l±g 


and 


cot(w-^-^-)  = 


T*±*J 


(25) 


As  a  consequence,  Eqs.(24)  reduce  Eqs.(A9)  and  (A  10)  to  the 
single  equation 


g+  sin0+  g_(Tosin0  Q 
1+g  1-g 


(26) 


As  usual,  the  vanishing  of  the  determinant  of  the  coefficient 
matrix  in  Eq.(17),  viz., 

Dy  (n*ff  )  =|M)'  (neff  )| =  0  .  (19) 

determines  the  effective  indices  rc^for  the  allowed  modes.  It  is 
instructive  to  note  that  in  the  limit  of  zero  chirality  the  matrix  My 
can  be  reduced  to  block  form  and  the  determinant 


Dy^neff\ 


M 


=  DQ(nejf )  reduces  to 


and  Eqs.(25)  reduce  Eqs.(Al  1)  and  (A  12)  to  the  single  equation 

D+  Sl sinb+d  -  ) ,  g-  g£ sinlw-rf  -£_ )  Q  (2?) 

1  +  h  \-h 

Requiring  the  determinant  of  the  coefficients  of  B+  and  B'  to 
vanish  in  Eqs.(26)-(27)  gives  another  condition  that  can  be 
written 


Do  (neff  )  “  II ^TM  *  (20) 


where 


MTe  - 


°o 

C 


-1 

-s 


Mm  - 


ro°o 

Cr 


-1 

-Sr 


(21) 


fl  +  h}_  +  Sq8 

,l-ftJ  SrJ-Sog’ 

(28) 

=<7q  sin(w±(i)-cos(w±^), 

+  +  +  + 

S/q"  =  q)(TQ  sin(w_rf)-cos(M“i/) , 


The  matrix  elements  in  Eq.(21)  are  devoid  of  ±  superscripts  and 
correspond  to  the  respective  quantities  defined  in  Eq.(18)  but 
with  y=0.  The  separate  vanishing  of  the  determinants  in  Eq.(20) 
gives  rise  to  the  well-known22  TE  and  TM  mode  equations,  which 
will  be  further  discussed  below. 

Eq.(19)  only  yields  the  allowed  effective  indices  n ^  To 
determine  the  remaining  parameters,  as  well  as  an  alternate 
approach  to  determining  n ^  we  return  to  Eqs.  (A9-A12)  and 
require  the  vanishing  of  the  determinant  of  the  coefficients  of  B+ 
and  B'  in  Eq.(A9)  taken  with  (A  10)  and  of  the  determinant  of  the 
coefficients  of  B+  and  B  in  Eq.(ll)  taken  with  (12).  After  some 
algebraic  manipulation,  we  can  write  these  conditions  in  the 
respective  forms 

1  t  11 
j  cot (0+)  ^  co 

GO 
and 


J  _  COt(u+d -0  +  )  j^COt(M  d-<t>  ) 

where  Aq  =  -j(l-q))  and  =  -j(l-rs) .  These  equations  are 

identically  satisfied  upon  introduction  of  new  parameters  g  and  h 
such  that 


(22) 


(23) 


where  we  have  used  trig  identities  and  Eqs.(24)  together 
withaQcr^  =o,^(To  ,  which  follows  from  their  definitions  in 

Eqs.(18).  Eqs.(24),(25)  and  (28)  represent  five  equations  for  the 
5  unknowns  t))*,  g ,  h ,  and  n $  and  contain  all  conditions  on  these 
parameters  imposed  by  the  boundary  conditions. 

Equations  for  the  subset  (g,  neff)  can  be  obtained  by 
eliminating  ty*  in  Eqs.(24)  and  (25)  and  also  solving  (28)  for  h  as 
a  function  of  g  to  arrive  at 


+  ,  -1 
u~d=  cot 


fa****] 

-if  +  r,±h^\ 
+  COt  1  0T  — - 

0  *±*  J 

1***1 

+  m~7i 


(29) 


fi(gtnej ff) 


(Sr£  ~Sr\)  )  +  (Sq  +5q  )g 
(SrQ  +SrJ)  +  (SQ  -S())g  ’ 


(30). 


where  m±  are  integers  that  are  not  necessarily  equal.  Using 
Eq.(30),  the  two  equations  in  (29)  can  be  simultaneously  solved23 
for  g  and  neS  after  which  the  phases  (J)*  can  then  be  found  from 
(24).  Eqs.(29)  reduce  to  the  equation24  for  the  TM  modes  of 
achiral  media  in  the  limits  (y-»0,  g-*0)  =»  /i — >0,  although  it  also 
has  solutions  with  TE  mode  limits  as  y-»0  as  explained  in  the 
following  text.  For  each  value  of  m  in  Eq.(29)  there  are  two 
solutions  for  neff:  one  with  |g|<land  one  with  |g|  >  1 .  The 

former  has  the  above-mentioned  TM  mode  limit  for  y-»0,  while 
the  latter  has  a  TE  mode  limit  for  y-»0.  Formally,  the  TE  mode 
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limit  of  Eqs.(29)  and  (30)  is  (y-^0,  g— >  ±°°)  =>  h^>  ±°°.  Because 
efficient  numerical  solution  of  Eq.(29)  requires  good  initial 
guesses  for  the  unknowns  n ^  and  g,  it  is  advantageous  for  the 
modes  with  a  TE  limit  ( viz .,  those  with  |g|  >  1 )  to  change  to  new 
parameters 

g=l/g,  h=\/hy  (31) 

in  order  to  limit  the  parameter  space  in  which  to  search  for  a 
solution.  With  these  substitutions,  Eqs.(29)  and  (30)  take  the 
forms 

J  f  T.+\ 


u±d=  cot  1 


a. 


+  ^g_±l 

0  *±i 


+  cot 


+  rsh±l 
— 


h  ±1 


+  m~7t 


T,~  ,  (Sr0  +sf>  )g+(So  ~So) 

k(g'neff)  =  —i7 - T - — 

(Sr0+-S/o)g+(5++S0) 


(32) 


(33) 


and  the  parameter  space  for  solutions  with  a  TE  mode  limit  is 
then  restricted  to  |g|  <  1 .  In  both  cases,  values  of  neff  are  limited 

to  max(rt0,  ns)<  neff<  min(n+,  n.)  if  both  the  RHC  and  LHC 
components  of  the  modes  [see  Eqs.  (14)-(16)]  are  to  be  totally 
internally  reflected  as  required  for  completely  guided  modes. 
Once  neff  and  g  (or  g  )  are  found,  the  phases  (j)*  can  be  obtained 

from  Eq.(24)  and,  selecting  B~  as  the  arbitrary  amplitude,  the 
remaining  constants  A*,  and  C4  can  be  obtained  from 
Eqs.(26),  (30),  (A1-A2),  and  (A5-A6): 


B+  =B' 


8+Jb 


8-' D  ) 


COS  0" 


COS0^ 


A~=B  cos0" 


g~> 0 


C±=B~cosiu~d-(j)  (/i±^ 

l  ^  J 


(35) 


(36) 


Note  that  for  modes  involving  solutions  g  of  Eq.(32),  the  phases 
(j)*,  and  constants  A*,  #+,  and  C4  can  be  found  from  Eqs. (24)  and 
(34-36)  simply  by  first  calculating  the  corresponding  value  of  g 
using  g  =  1/g  . 


parameters  (as  well  as  their  counterparts  g  and  h )  are  measures  of 
the  eccentricity  of  the  polarization  ellipse  for  the  transverse 
electric  field  in  the  top  and  bottom  claddings. 

From  Eqs.(3),  (8),  and  (9),  the  transverse  electric  field 
is  given  by 


Er  =  — 


ty  +  (y)+V~(y)]x  +ineff 


x  txp(-ik0neffz) 


v+(y)  v  (y) 


(37) 


which  is  recognized  to  be  elliptically  polarized  in  general,  with 
major  axes  lying  along  either  the  x-  or  y-axes.  Using  Eq.(35)  in 
(37),  the  ratio  of  the  y-component  to  the  jc-component  everywhere 
in  the  top  cladding  (achiral)  is  given  by 


'  _  •  neff  ^  _  .  neff  1 
|y>0  n0  g  ng  g  ‘ 


(38) 


Noting  that  nejf  jng  « 1 ,  especially  far  from  cutoff,  the 

polarization  state  of  the  transverse  field  in  the  top  cladding  can  be 
inferred  from  the  value  of  g  as  follows: 


>+l) 


TM: 

TE: 

g-»±“ 

(or 

(34) 

~RHC: 

(or 

-LHC: 

g- »-l 

(or 

(39) 


The  parameter  h  gives  the  polarization  state  of  the  transverse  field 
in  the  bottom  cladding  in  a  similar  fashion  because  the  use  of 
Eq.(36)  in  (37)  gives 


_  .  nejf  fis  _  ,  neff  1 
\y<-d  ns  h 


n8  h 


(40) 


In  the  core,  the  eccentricity  of  the  polarization  ellipse 
varies  as  a  function  of  y.  The  polarization  state  there  is  seen  to 
depend  on  and  as  follows:  ZT  =  0=>RHC,  #*  =  0=> 
LHC,  asp-»0]=»TE,  [B~-*-B+  as  y-»0]  =>  TM. 

From  Eqs. (24)  and  (26),  we  can  derive  the  equation 


4.  DISCUSSION 

Polarization  of  the  Transverse  Electric  Field 

An  essential  step  leading  to  the  form  of  the  solution  to 
the  chiral  core  asymmetric  planar  waveguide  given  in  the 
previous  section  is  recognizing  that  two  of  the  constraints  arising 
from  the  boundary  conditions  can  be  written  in  the  particular 
forms  given  in  Eqs.(22)  and  (23),  which  can  then  readily  be  seen 
to  be  identically  satisfied  by  introduction  of  the  parameters  g  and 
h  as  in  Eqs.(24)  and  (25).  We  will  now  show  that  these 


B  + 
B~ 


=  sign 


(g± 1) 
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y-axis  =  major  axis  of  polarization  ellipse 
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Figure  2.  Plots  of  the  eccentricity  parameter  g  versus  chirality 
for  the  first  six  elliptical  modes  in  a  planar  waveguide  with  a  3 
pm  thick  core  and  ng  =  1 .59,  n0  =  1 .00,  ns  -  1 .5 1 .  Top  plot: 
modes  that  are  RH  elliptical  for  y  >  0  (and  LH  elliptical  for  y  < 
0)  evolving  from  TM  modes  at  y  =  0.  Bottom  plot:  modes 
that  are  LH  elliptical  for  y  >  0  (and  RH  elliptical  for  y  <  0) 
evolving  from  TE  modes  at  y  =  0. 


which  relates  the  polarization  state  in  the  core  to  the  parameter  g . 
In  the  limit  of  zero  chirality,  in  which  case  <Tq  <7q  ,  the  TE 

and  TM  mode  limits  are  obtained  for  g— >  ( g  — >  0 )  and  g— » 

0,  respectively.  For  the  allowed  modes  in  a  chiral  core 
waveguide,  the  transverse  electric  field  is  elliptically  polarized. 
When  the  amount  of  chirality  is  small  (as  discussed  more  fully 
later),  the  modes  are  nearly  TE  (or  TM)  because  the  major  axis  is 
much  larger  than  the  minor  axis.  Modes  with  a  predominant 
RHC  polarization  or  LHC  polarization  are  realized  when  g^>  ±1, 
as  can  be  seen  from  Eq.(41).  There  is  a  maximum  on  the  right- 
hand  side  of  Eq.(41)  as  g->  1,  so  that  B+fB~  »1  resulting  in 
predominantly  RHC  polarization,  and  a  minimum  as  g— >  -1,  so 
that  B*IB~  «  1  resulting  in  predominantly  LHC  polarization. 


Polarization  Dependence  on  Chirality 

In  an  achiral  film,  as  is  well  known,  the  modes  are 
transverse  electric  (TE)  and  transverse  magnetic  (  TM).  An 
example  of  the  transition  with  increasing  chirality  from  TE/TM 
modes  in  the  achiral  case  to  predominantly  right-  and  left-handed 
circularly  polarized  modes  is  demonstrated  in  Figure  2  for  a 
planar  waveguide  with  a  3  pm  thick  core  and  ng  =  1.59,  n0-  1.00, 
ns=  1.51.  These  values  are  illustrative  of  a  chiral  polymer  on  a 
glass  substrate  at  a  wavelength  of  633nm.  The  parameter  g, 
which  is  proportional  to  the  eccentricity  of  the  polarization 


Figure  3.  Effective  index  neff  vs  thickness  d  curves  for  TE  and 
TM  modes  in  an  achiral  (y  =  0)  asymmetric  planar  waveguide 
with  ng  -  1.59,  n0=  LOO,  ns=  1.51. 


ellipse,  is  plotted  versus  the  chirality  parameter  y  for  the  first  six 
elliptical  modes.  The  top  plot  shows  the  modes  that  are  RH 
elliptical  for  y  >  0  (and  LH  elliptical  for  y  <  0)  evolving  from  TM 
modes  at  y  =  0,  and  the  bottom  plot  shows  the  modes  that  are  LH 
elliptical  for  y  >  0  (and  RH  elliptical  for  y  <  0)  evolving  from  TE 
modes  at  y=  0.  For  convenience,  the  curves  are  labeled 
according  to  the  handedness  of  the  mode  for  positive  y.  These 
curves  are  solutions  of  Eqs.(29)  and  (30)  for  the  case  |g|<l,  and 
Eqs.(32)  and  (33)  for  |g|>l  .  The  horizontal  dotted  lines  identify 
the  values  g  =  +1  and  g  =  -1  corresponding  to  predominantly 
RHC  and  LHC  polarization,  respectively,  of  the  transverse 
electric  field.  Note  that  for  a  given  chirality,  the  higher  the  order 
the  more  elliptical  (less  circular)  are  the  modes.  The  polarization 
eccentricity  for  a  given  elliptical  mode  also  depends  on  the  core 
thickness,  as  discussed  in  the  next  section. 

Thickness  dependence  and  cutoff  conditions 

Figures  3  and  4  show  the  mode  effective  index  as  a 
function  of  increasing  film  thickness  for  an  achiral  (y  =0)  and  a 
chiral  (y  =  4.5x10 5  pm)  waveguide,  respectively.  In  the  well- 
known  case  of  the  achiral  waveguide  (Figure  3),  note  that  the 
effective  indices  for  TE  and  TM  modes  of  a  given  order  get 
closer  together  as  the  thickness  increases  and  there  are  no 
crossover  points.  In  the  chiral  waveguide  (Figure  4),  on  the  other 
hand,  the  right-handed  elliptical  (RHE)  and  left-handed  elliptical 
(LHE)  modes  of  a  given  order  separate  as  the  core  gets  thicker 
and  the  RHE  modes  cross  over  LHE  modes  at  numerous 
thickness  values,  as  shown  more  clearly  in  the  inset  in  Figure  4. 
In  addition  to  the  usual  cutoff  at  a  minimum  thickness,  the  LHE 
modes  all  experience  a  second  cutoff  with  increasing  thickness 
determined  by  the  chirality  and  refractive  indices.  This  mode 
separation  and  mode  crossover  with  increasing  thickness  were 
also  found  in  the  derivations  in  Refs.  12  and  18,  but  neither  of 
these  previous  works  discussed  the  existence  of  the  second  cutoff 
for  completely  guided  LHE  modes. 

Just  like  the  achiral  case,  the  cutoff  at  minimum 
thickness  for  the  LHE  and  RHE  modes  occurs  when  neff  is  at  its 
minimum  value,  i.e.,  the  larger  of  no  and  ns.  Assuming  no<ns,  the 
cutoff  thickness  is  found  by  simultaneously  solving  for  the  cutoff 
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Figure  4.  Effective  index  ntS  vs  thickness  d  curves  for  RHE 
and  LHE  modes  in  a  chiral  (y  =  45  pm)  asymmetric  planar 
waveguide  with  ng  =  1.59 ,  /!#  =  1.00,  ns-  1.51.  The  inset  is  a 
blowup  emphasizing  the  mode  crossover  points  and  second 
cutoff  of  the  LHE  modes. 


Figure  5.  A  plot  of  the  effective  index  (left  axis)  and 
eccentricity  parameter  g  (right  axis)  for  the  RHEo  mode  versus 
thickness  illustrating  the  variation  in  eccentricity  in  the  vicinity 
of  the  intersections  of  the  RHEo  mode  effect  index  (left  axis) 
with  the  LHE  modes  (faint  lines)  for  the  same  waveguide 
parameters  as  in  Figure  4.  The  icons  included  to  help  visualize 
the  eccentricity  are:  vertical  double  arrow  -  TM,  horizontal 
double  arrow  -  TE,  ellipses  -  elliptical  polarization  with  x  or  y 
major  axis. 


values  gc  and  dc  with  n#  =  ns  substituted  in  the  two  equations 
(29).  (Because  arccot  («>)  =  0,  the  equations  do  not  contain  hc .) 
These  equations  can  be  conveniently  solved  in  two  steps.  First, 
gc  is  found  as  the  root  of 


Uc  (0r  +  Mk)  -  uc  (<t>c  +  Mn)  =  0  , 


(42) 


where 


"c  =*o^(rt±)2-rtj  • 


Pc  = cot  ' 


(l±s> 


|(n±)2-n|  /p±gc 


2  2 
-«o 


1  ±8c 


(43) 


and  M  represents  the  mode  number  M*  for  RHE  or  M~  for  LHE 
modes.  Once  gc  is  found,  the  cutoff  thickness  at  a  given 
wavelength  can  then  be  obtained  from 


Me  = 


,  +  w  + 

0  ~  +  M  ~7t 

^(n1)2  -n2 


(44) 


The  cutoff  value  hc  is  given  by 


RHE(Af=A/f)  =>  0<gc<l,  LHE(Af=AT)=>  -*°<gc<-ro 
y<0 : 

RHE(A/=AT)  =>  r0<  gc  <°°»  LHE(Af=AT)  =>  -1  <  gc  <0  . 

(46) 

For  those  cases  in  Eqs.(46)  where  \gc |>1,  it  is  advantageous 
numerically  to  solve  instead  for  gc  using  0/  in  place  of  0<r  in 

Eqs.(42)  and  (44),  where  0C*  has  the  same  definition  as  0^  in 
Eq.(43)  but  with  gc  replaced  by  Mgc . 

The  second  cutoff  at  maximum  thickness  for  LHE 
modes  (RHE  modes  when  y<0)  arises  because  both  of  the 
equations  (corresponding  to  the  two  signs)  in  Eq.(32)  must  be 
satisfied  simultaneously  for  a  given  mode  and  the  maximum  real 
value  of  nfff  for  completely  guided  modes  is  ntff  -  n+  (when  y  > 
0).25  No  matter  how  dominant  one  handedness  of  circularly 
polarized  waves  becomes  over  the  other  handedness,  both  types 
of  waves  must  be  present  in  order  to  satisfy  the  boundary 
conditions.  [Recall  that  both  LHE  and  RHE  waveguide  modes  are 
comprised  of  the  total  field  from  two  RHC-polarized  waves  and 
two  LHC-polarized  waves,  as  in  Eqs.(14)  and  (15),  and  are  the 
refractive  indices  associated  with  these  waves.]  Thus,  even  for 

LHE  modes  «+  must  be  real,  which  requires  that  nejf  <  . 


K  =  rs 


cos {ucdc)(r0  +  gc)  |  1 
cos  (u+dc)(r0-gc) 


(45) 


In  Eqs.  (42)-(45),  the  ranges  of  solutions  for  gc  are 
y>0 : 
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Figure  6.  Relationship  of  the  substrate  eccentricity  parameter  h 
(short-dashed  line)  to  the  upper  cladding  eccentricity 
parameter  g  (solid  line)  as  a  function  of  thickness  for  the  RHEq 
mode  of  Figures  4  and  5.  The  inset  shows  h~g  when  nf  -  m+ 
is  even  and  h  ~  1/g  when  nf  -  nf  is  odd.  Also  shown  for 
reference  is  the  core  amplitude  ratio  ZT/ZT  versus  thickness. 


Noting  that  the  cutoff  here  occurs  also  with  nearly 
circular  polarization,  we  put  n ^  =  n+  together  with  g  =h  =  -1  in 
Eq.(32)  to  obtain 


dc2(M~)  = 


cot 


-1 


+  cot 


d-<5). 


l(02-(n  +  )2  r0  +  l 


(n+)2  -n2 


( n  )2  ~(«+)2  rs  +1 
(n+)2-n2  2 


(47) 


+  M  n 


U0yj(n-)2-(n  +  )2 


Because  8  «  1 ,  an  approximate  value  for  this  cutoff  thickness  is 
given  by 


dc2(M ")* 


1 


2kQnR 


(M  +1)_ 

,  +  ,t)  j  1  +  ri'  1 

#  'T 

yj\-rs  [ 

.(48) 


From  this  expression  it  can  be  seen  that  in  the  achiral  limit 
(whence  8  0),  this  second  cutoff  thickness  becomes  infinite, 

consistent  with  the  nonexistence  of  such  a  second  cutoff  in  the 
achiral  case.  The  maximum  thickness  cutoff  of  the  LHEq  mode 
for  the  example  in  Figure  4,  which  has  y=4-5xl0'5  pm,  occurs  at 
slightly  less  than  3  pm.  This  value  of  y  corresponds  to  a  bulk 
rotatory  power  of  642  deg/mm  at  a  wavelength  of  633nm.  For 
reference,  a  bulk  rotatory  power  of  30  deg/mm  (y=2.lxl0~6  pm) 


d  (^im) 


Figure  7.  Detail  of  the  crossover  of  mode  RHEq  (dotted  line) 
with  LHEj  (solid  line)  in  Figures  4-6.  Top  plot:  effective  indices 
versus  thickness.  Bottom  plot:  eccentricity  parameter  g  for  RH 
polarization  {Mg  for  LH  polarization)  versus  thickness.  In  the 
narrow  transition  region,  indicated  by  vertical  dotted  lines,  only 
LH  polarization  exists. 


results  in  a  maximum  thickness  cutoff  of  16.7  pm  for  the  LHEq 
mode. 

At  this  point,  we  note  the  connection  between  our 
solutions  and  the  approximate  solutions  for  the  asymmetric  planar 
waveguide  given  in  Ref.  18.  We  can  obtain  the  RHC  mode 
equations  of  their  work  by  putting  g= 1,  h=  1,  rs  =1,  and  solve 
Eq.(29)  with  the  plus  sign  while  ignoring  the  corresponding 
equation  with  a  minus  sign.  Similarly,  their  LHC  modes  can  be 
obtained  by  setting  -1,  /i=  -1,  rs  =1  and  solve  Eq.(29)  with 
the  minus  sign  while  ignoring  the  companion  equation  with  a  plus 
sign.  This  gives  reasonably  good  values  for  the  effective  indices 
of  the  RH-modes  in  general  and  for  the  LH  modes  between  the 
low  and  high  thickness  cutoff  values.  However,  the 
approximations  used  in  Ref.  18  resulted  in  missing  the  second 
thickness  cutoff  for  the  LHE  modes  (y>0)  and  ignoring  the  fact 
that  higher  order  modes  have  differing  degrees  of  elliptical 
polarization  (see  Figure  2). 

Eccentricity  and  mode  crossover  points 

In  this  section  we  examine  the  behavior  of  mode 
eccentricity  with  increasing  thickness,  the  relationship  between 
the  integers  m*  and  the  mode  order  numbers  AZ*  ,  and  the  details 
of  the  mode  crossover  points.  In  Figure  5  is  plotted  the  effective 
index  (left  axis)  for  the  RHEq  mode  along  with  the  eccentricity 
parameter  g  (right  axis)  as  a  function  of  thickness.  With 
increasing  thickness  d ,  the  effective  index  neff  for  RHEq  crosses 
over  successively  higher  LHE  modes,  viz.,  LHEIt  LHE2,  etc.  To 
obtain  this  curve  for  mode  M ¥ =0,  Eqs.(29)  and  (30)  are  solved, 
starting  from  cutoff,  with  nf  =  0  and,  initially,  nf  =  0.  The 
eccentricity  parameter  g  for  the  upper  cladding  starts  out  small  in 
this  example  (nearly  TM,  with  the  y-axis  the  major  axis  of  the 
polarization  ellipse)  but  increases  to  g  =  1  (circular  polarization) 
as  d  increases.  There  is  then  a  n  discontinuity  in  arccot  and  nf 
increases  by  1  to  nf  =  1.  With  further  increase  in  thickness,  the 
;t~axis  becomes  the  major  polarization  axis  (because  g  >1)  and 
eventually  g  (TE)  as  RHEq  crosses  over  LHEj.  After 
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Figure  8.  Effective  index  neff  vs.  thickness  k$  mode  curves  for  a 
2  fim  thick  core  with  ng  =  1.59,  n0=  1.00,  and  ns  -  1.57.  The 
LHE  modes  are  cut  off  by  the  line  neff  =  n+(ko ),  while  the  RHE 
modes  are  asymptotic  to  this  line  undergoing  a  high  frequency 
cutoff  when  n+(k o)  — >  ns . 

crossing  over  LHElt  the  sequence  repeats,  starting  with  g  =  0  ( 
TM),  then  RH  elliptical  with  major  y-axis,  increasing  to  g  =  1 
(circular,  where  nf  is  bumped  up  to  nf  =  2),  then  RH  elliptical 
with  major  jc-axis,  and  finally  g  (TE)  as  the  LHE2  crossover 
approaches.  Polarization  ellipse  icons  are  included  in  Figure  5  to 
help  visualize  this  sequence.  In  general,  for  an  RHE  mode  Af\ 
we  fix  nf  -  M*  in  Eqs.(29)  and  (30)  [or  Eqs.(32)  and  (33)]  and  nf 
varies  according  to  nf=  M \  Nf  +  1,  Nf  +  2, ...  with  increasing  d. 

The  polarization  parameter  h,  which  is  proportional  to 
the  eccentricity  in  the  substrate,  satisfies  h~g  when  nf  -  nf  is 
an  even  integer,  and  h  ~  Mg  when  nf  -  nf  is  an  odd  integer. 
This  is  shown  in  Figure  6,  where  h,  g,  nf,  and  B+/B~,  are  plotted 
versus  thickness  for  the  RHEq  mode.  The  inset  shows  h  plotted 
directly  against  g  . 

A  similar  sequence  of  events  occurs  for  the  LHE  modes 
with  the  differences  that  now  g  is  negative  and  a  given  LHE 
mode  only  crosses  over  RHE  modes  of  lower  mode  number  as  the 
thickness  increases.  In  general,  for  LHE  mode  NT,  we  fix  nf  -M 
and  nf  varies  according  to  nf  =  AT,  M~  -  1,  NT  -  2, ... ,  -1 
with  increasing  thickness,  with  cutoff  occurring  while  nf  =  -1  . 

It  would  seem  that  there  is  an  abrupt  polarization 
discontinuity  at  the  cross-over  points  in  Figure  5.  However,  on 
closer  examination,  it  turns  out  that  there  is  a  transition  region  at 
the  mode  crossovers,  as  illustrated  in  Figure  7  with  vertical  dotted 
lines  for  the  RHEo/LHEj  crossover.  Inside  this  small  region,  both 
modes  are  LH  elliptical  -  one  with  an  jc-major  axis  and  one  with  a 
y- major  axis.  The  effective  index  curve  for  RHEq  merges 
continuously  from  the  left  into  the  LHEj  curve,  and  the 
eccentricity  for  RHEq  varies  continuously  from  RH  elliptical, 
with  jt-major  axis,  to  TE  (l/g=0),  to  LH  elliptical  with  jc-major 
axis.  The  effective  index  curve  for  LHE!  merges  continuously 
from  the  left  into  the  RHEo  curve,  and  the  eccentricity  for  LHEj 
varies  continuously  from  LH  elliptical,  with  y-major  axis,  to  TM 
(g=0),  to  RH  elliptical  with  y-major  axis.  Although  this  behavior 


Figure  9.  Eccentricity  parameter  g  for  the  first  three  RHE  modes 
(top  plot)  and  Mg  for  the  first  three  LHE  modes  (bottom  plot) 
versus  thickness  for  chirality  y=045  picometers  (solid  line)  and 
Y=4.5  picometers  (dotted  line)  in  a  symmetric  waveguide  with 
ng  =  1.59  and  n0=ns-  1.57.  The  values  at  cutoff  [see  Eq.(52)] 
for  all  but  the  zeroth  order  modes  are  gc  =  ±rs,  where,  in  this 
example,  rs  =  0.975  ;  note  the  initial  rapid  falloff  and  subsequent 
rise  for  d>dc . 

contradicts,  somewhat,  the  designation  RHE  inside  this  region, 
we  maintain  the  RHE/LHE  mode  designations  because  the 
crossover  transition  regions  are  very  small. 

In  general,  the  thickness  values  for  the  transition  region 
edges  for  the  RHE^  mode  at  the  LHE^-  crossover  region  can  be 
obtained  from  Eqs.(29)  using 

AT  -  Mf  -  odd  integer 

low-d  side:  1/g  =  0,  h  =  0 
high-d  side:  g  =  0,  Mh  =  0; 

(49) 

AT  -  A/*  =  even  integer 

low -d  side:  Mg  =  0,  Mh  -  0 
high-dside:  g  =  0,h  =  0. 

These  substitutions  in  Eq.(29),  together  with  nf  -M~  and  nf  = 
M+,  yield  two  equations  with  unknowns  neg and  d.  They  can  be 
straightforwardly  solved  by  first  eliminating  d,  solving  the 
resulting  equation  numerically  for  neff,  and  then  using  either  of  the 
original  equations  to  find  d. 

Frequency  dependence 

Ref.  18  predicted  two  frequency  cutoff  values  for  the 
right-handed  modes.26  We  will  show  that  both  RHE  and  LHE 
modes  have  a  low-frequency  and  a  high-frequency  cutoff  for 
completely  guided  modes,  although  for  different  reasons.  The  low 
frequency  cutoff  is  the  usual  one  that  corresponds  to  the 
maximum  wavelength  Xc  for  which  a  waveguide  of  given 
thickness  d  can  still  support  the  mode  A/*  of  interest.  This  can  be 
found  by  simultaneously  solving  for  the  cutoff  values 
k Qr  =  2n/Xc  and  gc  in  the  two  equations 
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uid  =0;r  +m±n , 


(50) 


where  u *  and  0*  are  given  by  Eqs.(43),  taking  into  account  that 
&=kongy  and  n±  ~  ng /(l±5)are  wavelength  dependent.  At  this 
cutoff  we  put  m+~m~=M+  for  RHE  modes  and 

m+  =  m~  -M~  for LHE modes. 

Starting  from  these  cutoff  values,  Figure  8  shows  the 
effective  index  for  the  first  two  LHE  and  RHE  modes  plotted  as  a 
function  of  Icq  =  2 7ijX  for  a  2  pm  thick  core  with  ng  =  1.59,  n0- 
1.00,  and  ns=  1.57.  In  the  interest  of  simplicity,  the  frequency 
dependence  of  the  chirality  parameter  y  and  of  the  average  index 
ng  is  not  taken  into  account  in  this  plot.  For  reference,  the 
frequency  dependence  of  the  effective  index  for  the  first  two  TE 
modes  are  plotted  for  the  achiral  guide  with  a  core  index  of  ng. 
(The  TM  modes  are  virtually  indistinguishable  from  the  TE 
modes  on  this  scale.)  Just  as  there  is  a  maximum  thickness 
cutoff  for  completely  guided  LHE  modes,  there  is  a  high 
frequency  cutoff  that  occurs  as  —>  n+(ko)  for  these  modes.  The 
numerical  value  can  be  obtained  by  finding  the  root  koc2  that 
satisfies  Eq.(48)  [after  replacing  dc2  with  the  desired  thickness  d\ 
together  with  (6).  The  corresponding  left-handed  polarized  modes 
in  Ref.  18  have  no  high  frequency  cutoff  and  continue 
indefinitely  with  increasing  fa 

A  high  frequency  cutoff  for  the  right-handed  modes18 
seen  in  Figure  8  arises  because  in  order  for  completely  guided 

modes  to  exist  the  effective  index  is  limited  to  ns  <  nejf  <  n+ 

(here  assuming y>0  and  n5>n0).  As  n+  =ng /(I  +  5)  decreases 
with  increasing  5  (decreasing  X),  there  will  be  some  sufficiently 
short  wavelength  such  that  the  limit  n+  ->  ns  =»  8  *  (1  -  rs)j 2 

is  reached,  and  no  guided  modes  are  allowed. 

The  eccentricity  parameter  g  also  undergoes  dramatic 
changes  as  a  function  of  frequency  at  a  fixed  thickness,  like  those 
changes  shown  in  Figs.  (5)  and  (6)  for  the  case  of  changing 
thickness  and  fixed  frequency. 


The  symmetric  chiral-core  planar  waveguide 

As  the  refractive  index  rt0  of  the  upper  cladding 
approaches  n5  for  the  lower  cladding,  the  ranges  of  thickness 
values  where  m~  -  m+  is  an  odd  integer  (for  example,  the  regions 
in  Figure  6  where  h~  Mg)  become  narrower  and  narrower  until, 
in  the  limit  n0  =  ns>  only  even  integer  differences  exist  and  h  -  g 
everywhere.  In  this  limit,  Eq.(29)  reduces  to 

u±d  -ifa  +m±7i ,  (51) 

with  (J)*  given  by  Eq.(24).  [Note  that  Eq.(25)  is  now  equivalent  to 
(51)].  When  Eq.(51)  is  solved  for  <)>*  and  substituted  in  Eq.(l  1), 
the  solutions  to  Eqs.(10)  in  the  core  take  the  forms 

Vcore(y)=B±  coslp±(y  +  d/2)+  m*  71/2),  0<y<-d,(52) 

which  are  easily  seen  to  separate  into  even  or  odd  functions  (viz., 
pure  cosines  or  sines)  depending  on  whether  m±  are  even  or  odd 


integers,  respectively.  As  the  waveguide  thickness  increases  for 
RHE  modes  (m+  =  M  +),  the  integer  nf  progresses  through  the 
series  m"=  A/*,  At  +  2,  Af*  +  4,  ...,  with  the  change  to  the  next 
higher  value  occurring  as  g  passes  through  +1.  For  the  LHE 
modes  (m~  =  M  “  ),  m+  varies  according  to  tn  =  Af  “  -  2  (for 
M~  *  0),  M~t  M~  -  2,  M~  -  4, ...  with  the  change  to  the  next 
value  occurring  as  g  passes  through  -1;  the  second  cutoff  at 
maximum  thickness  occurs  at  m+  =  0  for  even  modes  and  m+  =  -1 
for  odd  modes. 

The  cutoff  at  minimum  thickness  for  the  symmetric 
waveguide,  given  by  Eq.(44)  with  0^  —■> 0  as  — >  ns,  is  the 
same  as  found  in  Ref.  12.  As  in  the  achiral  case,  the  zeroth  order 
modes  persist  down  to  d= 0.  To  obtain  the  values  of  g  at  cutoff, 
we  use  Eq.(51)  with  (43)  in  the  limit  ns  =  n0.  For  the  cutoff 

of  RHE  modes  (g> 0,  we  have  0^  ->  0  and  fa  ±  0  , 

while  0c  *  0  and  0^  -» 0  for  the  cutoff  of  LHE  modes  (g< 0, 

For  5*0  we  get 


where 


Sc  = 


-G±^G2  +  rs,M±  =0 
±rs  ,  M±*0 


G  = 


*1 

5(1-A,)  ' 


(53) 


(54) 


It  is  notable  that  for  Af*  *  0,  the  cutoff  value  of  the  eccentricity 
parameter  g  is  independent  of  the  chirality  parameter  y,  so  long  as 
y*0.  This  is  illustrated  in  Figure  9,  which  shows  the  behavior  of 
g  near  cutoff  for  the  first  three  RHE  and  LHE  modes.  In  addition, 
as  follows  from  Eq.(41),  ZT= 0  at  an  RHE  cutoff  and  B+= 0  at  an 
LHE  cutoff.  The  second  cutoff  at  maximum  thickness  for  the 
LHE  modes  in  the  symmetric  case  is  given  by  Eqs.(47)  and  (48) 
with  r0  =  r5. 

5.  CONCLUSION 

In  this  paper  we  have  presented  the  general  solution  in  Eqs.(29) 
and  (30)  for  modes  in  an  asymmetric  planar  waveguide  with  an 
isotropic  chiral  core  in  terms  of  a  pair  of  parameters  g  and  h 
related  to  the  eccentricity  of  the  polarization  ellipse  for  the 
transverse  electric  field  [see  Eqs.(38)  and  (40)].  The  symmetric 
waveguide  is  discussed  as  a  limiting  case.  This  formulation 
provides  insight  into  the  transition,  with  increasing  chirality  of 
the  core,  from  TE/TM  modes  in  the  achiral  case  to  modes  with 
nearly  right-handed  and  left-handed  circular  polarization.  The 
polarization  of  the  transverse  electric  field  for  the  modes  is 
elliptical  in  general  with  RHE  (right-handed  elliptical)  modes 
evolving  from  TM  modes  and  LHE  modes  evolving  from  TE 
modes  as  the  chirality  parameter  y  increases  from  zero.  The 
eccentricity  of  the  polarization  of  the  transverse  electric  field 
varies  dramatically  with  thickness  and  with  frequency. 
Expressions  for  the  minimum  thickness  cutoff  values  are  given  in 
Eqs.(42)-(44)  for  the  asymmetric  waveguide  and  in  Eqs.(44)  and 
(53)  for  the  symmetric  waveguide.  For  materials  with  positive 
y,  completely  guided  LHE  modes  (RHE  modes  for  negative  y) 
experience  a  second  thickness  cutoff,  given  by  Eq.(48),  so  that 
total  internal  reflection  of  both  the  RHC  and  LHC  components  of 


9 


these  left-handed  modes  occurs  only  between  a  minimum  and 
maximum  thickness.  Correspondingly,  there  are  also  minimum 
and  maximum  frequency  cutoffs.  Although  the  existence  of  a 
maximum  frequency  cutoff  for  the  RHE  modes  has  been  reported 
earlier18,  the  second  cutoff  conditions  have  heretofore  not  been 
reported  for  the  LHE  modes.  The  unique  polarization  properties 
of  chiral  waveguides  discussed  in  this  paper  may  lead  to  new  and 
useful  device  designs. 

Appendix:  Boundary  Conditions 

The  continuity  of  the  tangential  components  of  E  and  H  ,  viz., 
Ex  ,  EZ,HX, and  H z  at  the  two  interfaces, y  =  0  and y--d , 
gives  rise  to  a  set  of  8  equations  -  one  less  in  number  than  the 
nine  unknowns  A±y  B±y  C4,  0*,  and  n eS-  leaving  arbitrary  the 
overall  amplitude  of  the  lightwave.  From  Eqs.(3)  together  with 
(8),  (9),  and  (1 1),  one  gets 

E1  at  v =0  : 

B+cos<p+  +B~cos<t>-  =A+  +A~  (Al) 

H1  at  v=0 : 

#+cos0+-tf“cos0”  -^{a+  -A~)  (A2) 


while  elimination  of  C4  in  the  last  four  equations  gives 


B+ 


|cx+  sin(u+d  -0  +  )-cos(w+rf  -0+)j 
+  B~  [rj  sin {u~ d  ~<p~)~  cos (u~d  -0_  )]=  0 


(All) 


and 


|r5cr4  sin(w+d -0  +  )-cos(w+d -0  +  )j 
-  B~ [rvcx~  sin(u~~d  -0_) -cos(u~d  -0_)]  =  0 


(A  12) 


As  discussed  in  the  body  of  this  paper,  these  last  four  equations 
determine  the  effective  indices  neff ,  the  phases  0*,  and  the  ratio 

b+/b~  . 
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Ez  at  v=0 ; 

B+fa+sintUB^crjsinr^-A-  (A3) 

iL at  v~fl : 

sin0+  +B~<Jq  sin0”  =A+  +  A~  (A4) 

Ex  at  y--d : 

B+ cos(u+d-<j)+)  +  B' cos(u~d-<j)~)  =  C+  +C“  (A5) 

H1  at  y=~d : 

B+  cos (u+d  -(p+)-B~  cos (u~d  -0_)  =  ^7(C+  -C”) 

(A6) 

Ez  at  v  ~-d  : 

B+ sin(w+^/  -(p  +  )-  B~ sin(w_af  -0”)  =  C+  -C~ 

(A7) 

Hz  at  \=-d\ 

B+o+  sin (u+d  -0+)+  B~aJ  sin(w_£/  -0_)  =  C+  +  C~ , 

(A8) 

where  cTq  ,  erf ,  /ft ,  and  rs  are  defined  in  Eq.(  1 8). 

Eliminating  A±  in  the  first  four  equations  gives 

B+(i Tq  sin0+  -cos0  +  )  +  fi_(ao  sin0“  -cos0“)  =  0  (A9) 

and 

S+(r0<T0+  sin0+  -cos0  +  )-fi“(q)tTo  sin0~  -cos 0“)  =  0  , 

(A10) 


REFERENCES 

1  D.L.  Jaggard,  A.R.  Mickelson,  and  C.H.  Papas,  “On 
electromagnetic  waves  in  chiral  media,”  J.  Appl  Phys.  19, 
211(1979). 

2  A.  Lahtakia,  V.K.  Varadan,  and  V.V.  Varadan,  Time-Harmonic 
Electromagnetic  Fields  in  Chiral  Media ,  Lecture  Notes  in 
Physics  Series  335  (Springer- Verlag,  Berlin,  1989). 

3  N.  Engheta  and  D.L.  Jaggard,  “Electromagnetic  chirality  and  its 
applications,”  IEEE  Trans .  Antennas  and  Propagation 
30,6(1988). 

4  V.K.  Varadan,  A.  Lakhtakia,  and  V.V.  Varadan,  “Propagation  in 
a  parallel  plate  waveguide  wholly  filled  with  a  chiral  medium,”  7. 
Wave  Mater.  Interact  3,  267(1988). 

5  N.  Engheta  and  P.  Pelet,  "Modes  in  Chiro waveguides,"  Optics 
Letters  14(11),  593(1989). 

6  P.  Pelet  and  N.  Engheta,  "The  Theory  of  Chirowaveguides,  ” 
IEEE  Trans.  Antennas  and  Propagation  38(1),  90(1990) 

7  L.  Zhang,  Y.  Jiao,  and  C.  Liang,  “The  dominant  mode  on 
parallel-plate  chirowaveguides,”  IEEE  Trans .  Microwave  Theory 
and  Techniques  42,  2009(1994). 

8  See  for  example,  A.  Yariv  and  P.  Yeh,  Optical  Waves  in 
Crystals  (New  York,  J.  Wiley  and  Sons,  1984),  section  4.9. 

9  For  a  recent  review,  see  M.M.  Green,  J.-W.  Park,  T.  Sato,  S. 
Lifson,  R.L.B.  Selinger,  J.V.  Selinger,  Angew.  Chem.  Int.  Ed.  38, 
3138  (1999). 

10  Highly  chiral  oligimers  have  been  realized,  for  example,  in 
metallocenes  [T.J.  Katz,  A.  Sadhakar,  M.F.  Teasley,  A.M. 

Gilbert,  W.E.  Geiger,  M.  P.  Robben,  M.  Wuensch,  and  M.D. 
Ward,  J.  Am.  Chem.  Soc.  115,  3182(1993).]  but  additional 
requirements  of  low  optical  loss  and  incorporation  into  polymers 
with  good  film-forming  properties  must  be  met. 

11  M.  Chien,  Y.  Kim,  and  H.  Grebel,  “Mode  conversion  in 
optically  active  and  isotropic  waveguides,”  Optics  Letters  14, 
826(1989). 


10 


12  M.  Oksanen,  P.K.  Loivisto,  and  L.V.  Lindell,  "Dispersion 
curves  and  fields  for  a  chiral  slab  waveguide, "  IEE  Proc.-H 
138(4),  327(1991). 

13  J.  Xiao,  K.  Zhang,  and  L.  Gong,  "Field  Analysis  of  a  General 
Chiral  Planar  Waveguide,"  Int.  J .  of  Infrared  and  Millimeter 
Waves  18(4),  939(1997) 

14  S.F.  Mahmoud,  “Guided  Modes  on  Open  Chirowaveguides,” 
IEEE  Trans.  Microwave  Theory  Techno L,  MTT-43,  205  (1995). 

15  K.M.  Flood  and  D.L.  Jaggard,  “Single-mode  operation  in 
symmetric  planar  waveguides  using  isotropic  chiral  media,” 
Optics  Letters  21,  474(1996) 

16  S.  Yamamoto,  Y.  Koyamada,  T.  Makimoto,  “Normal-mode 
analysis  of  anisotropic  and  gyrotropic  thin-film  devices  for 
integrated  optics,”  J.  Appl.  Phys.  43,  5090(1972). 

17  P.Pelet  and  N.  Engheta,  “Coupled-mode  theory  for  chiral 
waveguides,”  J.  Appl.  Phys .  67,2742(1990). 

18  S.V.  Demidov,  K.V.  Kushnarev,  and  V.V.  Shevchenko, 
“Dispersion  Properties  of  the  Modes  of  Chiral  Planar  Optical 
Waveguides,”  J.  Communications  Technology  and  Electronics 
44,827(1999).  Translated  from  Radiotekhnika  i  Elektronika  44, 
885(1999). 

19  For  a  discussion  of  various  forms  for  chiral  constitutive 
equations,  see  Ref.  2,  chap.  3. 

20  C.F.  Bohren,  “Light  scattering  by  an  optically  active  sphere,” 
Chem.  Phys.  Lett.  29, 458(1974)  and  “Scattering  of 
electromagnetic  waves  by  an  optically  active  spherical  shell,”  J. 
Chem.  Phys.  62,  1566(1975).  See  also  the  discussion  in  Ref.  2, 
chap.  7  and  Ref.  18. 

21  C.R.  Paiva,  A.L.  Topa,  and  A.M.  Barbosa,  “Semileaky  waves 
in  dielectric  chirowaveguides,”  Optics  Letters  17,  1670(1992). 

22  See  for  example,  P.K.  Cheo,  Fiber  Optics  and  Optoelectronics, 
2nd  Edition,  (Prentice  Hall,  New  Jersey,  1990),  chapter  3. 

23  Solutions  for  all  examples  in  this  paper  were  obtained  using 
built-in  routines  in  Mathcad  2000  Professional  (MathSoft,  Inc., 
Cambridge,  MA). 

24  The  TE  and  TM  mode  equations  are  normally  written  using 
arctangent  rather  than  arccotangent.  However,  the  arguments  of 
the  arctangent  are  positive  definite  and  arctan(jc)  =  arccot(//jc), 
for  *>0  as  is  true  for  the  achiral  case.  There  is  an  advantage  to 
using  arccotangent  in  the  chiral  case  because  singularities  of  the 
argument  cause  n  discontinuities  of  arccotangent  only  at  the 
endpoints  of  the  range  |g|<l,  while  if  arctangent  is  used,  there 
are  n  discontinuities  at  g  =  ±  r0  inside  this  range. 

25  For  a  discussion  of  semileaky  waves  in  chiral  waveguides,  see 
Ref.  21. 

26  The  modal  dispersion  equations  in  Ref.  18  are  written  using 
normalized  waveguide  parameters.  Here  we  have  refrained  from 
using  normalized  parameters  because  neither  the  normalization 
used  in  Ref.  18  nor  that  used  in  Ref.  12  can  simultaneously  be  used 
to  investigate  both  thickness  and  frequency  dependence. 


